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ABSTRACT 


The Naval Postgraduate School simulation model, QDYN&5, 
was used to investigate the sputtering of nitrogen reacted 
tungsten and molybdenum targets, where the nitrogen was 
placed in equilibrium positions both above and below the 
target surface plane. The cascades were initiated by 
bombarding the peace with argon ions with energies varying 
from 1 kev to 3 kev. The simulation results, when compared 
to experimental results obtained by Winters, indicate that 
the nitrogen position is not below the surface plane for 
either target but is located near the point where the adatom 
would be equidistant from its five nearest neighbors. 
Furthermore, the results show no multimer formation as a 
result of lattice fragmentation but do show dimer formation 
as a result Ot recombination above the target of 


individually sputtered atoms. 
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IT. OVERVIEW 


A. DEFINITION AND EARLY EXPERIMENTAL WORK 

Sputtering is the ejection of surface atoms from an 
object due to energy transfer from an energetic projectile 
ion to the atoms of the object. During each collision 
between the energetic ion and a target atom, some energy and 
momentum is given ae by the ion. If this transferred energy 
exceeds the atoms’ lattice binding energy, the struck atom 
(recoil atom) will be displaced out of its lattice site and 
will collide with other target atoms initiating a cascade. 
If a recoil atom is located at the surface of the ae 
and its acquired energy exceeds the surface binding energy, 
and the direction of motion is away from the target, that 
atom will be ejected. This cascade sequence continues until 
the initial energy ree been dissipated so far that further 
eoiliasions could not possibly result in target atom 
displacements. 

The sputtering process was first reported in 1852 by 
GROVE (Ref 1) when he noticed a deterioration of cathodes 
and a blackening of the glass of glow discharge tubes. Over 
the next several decades, a considerable amount of work was 
done in the field of sputtering. However, a lack of 
understanding of how various parameters such as ambient 


pressure, target purity and angle of incidence of the ion 


beam would affect the sputtering yield ,Y, ( defined as the 


number of ejected particles per incident ion) resulted in 
non-reproducible experiments. This effectively negated these 
early etforts. 

By 1902, GOLDSTEIN (Ret 2) presented evidence that 
sputtering was caused by positive ions impacting on the 
metal cathode. Since this was considered a deterioration and 
thus an undesireable effect at the time, work in the area 
was centered on Eon oleae the sputtering effect. However, 
some early researchers such as PLUCKER (Ref 3), recognized 
the possibilities for other experimental work of the effect 
in thin film manufacturing and removal. PLUCKER and other 
pioneers in sputtering laid the foundation which established 


sputtering as a field of quantitative experimental research. 


B. THEORETICAL AND EXPERIMENTAL DEVELOPMENTS 

The dependence of the sputtering yield on the ambient 
pressure was first demonstrated by PENNING and MOUBIS (Ret 
4) an 1940. They realized that the number of collisions 
between atoms ejected from the target and atoms of the 
ambient atmosphere of the vacuum system increased when the 
pressure was increased. They showed that this was causing 
backscattering of sputtered atoms into the target. To 
overcome this effect, they kept the ambient pressure no 
higher than 10-5 Torr. At this level, reproducible results 
were obtained for ion energies greater than 500 ev. For ion 


energies less than 500 ev, the sputtering action was 


iMsusteclenit to Maintain’ a ~clean surface. 

STARK (Ref 5 ), pursuing sputtering as an individual, 
atomic scale event, developed two separate models ot the 
process. His hot spot model was based on the assumption that 
the evaporation of target material from a very small area 
was caused by the high temperature created by the ion’s 
impact. VON HIPPLE (Ref 6) further developed this approach 
and expended a Ponecldoe ee amount ot effort attempting to 
formulate a sputtering theory based a this premise, stating 
that this approach was the only feasible avenue to address 
the complex statistics of the collision processes occuring 
in a sputtering event. It should be noted that although this 
approach has aeoeray been reexamined by KELLY and coworkers 
(Ref 7), it has not enjoyed the attention and development 
accorded to Stark’s other model, the “collision” model. 
In 1921, KINGDON and LANGMUIR (Ref 8) produced a_ successful 
application of Stark’s collision theory to the analysis of 
i0n induced desorption of monolayers. 

In the 1950’s, KEYWELL (Ref 9,10) attempted to formulate 
Stark’s collision model in terms of neutron transport 
theory. This and subsequent calculations by HARRISON (Ref 
11) were significant because probability concepts, in the 
form of collision cross sections, were introduced into the 
tneory™ Although important steps, these theories were 
not useful because they depended on a large number of 


unknown parameters; so quantitative yield values could not 


be calculated. In 1956, WEABNER (Ret 92) demonspeaceaeetne 
effect that the target’s crystalline structure could have on 
the sputtering process when he observed spot patterns while 
bombarding various monocrystals. These spots were atom 
ejection patterns which corresponded to the close-packed 
directions of the monocrystal target. 

Futher theoretical work by SILSBEE (Ref 13) led to the 
suggestion that eee, was transported through a lattice 
structure by momentum transfer along close packed 
directions. This phenomenon is called focusing and the 
transport mechanism is called a focuson. Extensions of 
Silsbee’s focusing concept were developed by THOMPSON (Ret 
14 ) and by experimental work by NELSON and THOMPSON (Ref 
15). While the focuson mechanism exists, subsequent 
investigations have shown that its contraibut lone to 
sputtering is small and that, by itself, it could not 
explain the spot patterns observed by WEHNER. By analogy, 
ROBINSON and OEN (Ref 16) expected that a channelling effect 
on an ion through the lattice would be found. This 
channelling effect would, of course, be dependent on the 
angle of incidence and the crystalline structure. They 
effectivly discovered this while studying ion penetration of 
a crystalline structure using a simulation model. This 
phenomenon was confirmed later experimentally by DAVIES et 
al. (Ref 17) 


Alternative theoretical approaches to the sputtering 
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problem culminated in the theories of SIGMUND (Ref 18) and 
THOMPSON (Ref 14) which consider ensemble averages over 
distribution functions. This approach suppresses any use of 
a target lattice structure in favor of an amorphous target. 
At higher energies, reasonably correct values of ite tre 
sputtering yield, the correct form for the ejected atom 
energy distribution function and the ejected atom angular 
Gast ri bution Menction from polycrystalline targets are 
obtained. This theory is relativly inaccurate at low ion 
energies since it is in this region that lattice effects are 
most pronounced and are of significant consequence. While 
refinements’) and Mae ouenere followed, attempts to reduce a 
many - bodied, multiple interaction event aco an 
analytically feasible model have still not been able to 


fully explain observed behavior. 


C. COMPUTER SIMULATION 

With the advent of high speed, sophisticated computers, 
iG became possible to simulate a sequence otf atomic 
collision events. The development of the simulation tool 
proceeded along two paths. Since the development of 
sputtering as a series of independent binary collisions was 
a natural outgrowth of STARK’S initial theories, this type 
of simulation was developed first. One of the principle 
philosophical drawbacks to this approach is that a complete 


understanding otf the physical model is assumed and that the 


i 


assumption that each event 15 separate from another can be 
justified. While this concept may be practical for higher 
initial energies, it is suspect for lower energies where 
multiple collisions are important. 

Starting in 1960, a different approach to the sputtering 
Simulation developed. It was recognized that, particularly 
at low energies, a cascade could not be treated aS a series 
of independent an Instead, a multiple interaction (MI) 
procedure was required. In this type of procedure, during 
a specific short period of time, many atoms of the target 
will simultaneously affect the cascade’s development. This 
type of Set was introduced by GIBSON, GOLAND, MILGRAM 
and VINEYARD (Ref 19) when they developed a computer 
Simulation model of the motions ot copper atoms ina 
crystalline target after one of the atoms had been struck by 
an incident ion. 

This time-step approach to sputtering simulation has, 
over the years, illuminated important aspects of the process 
which were masked by the ensemble averaging theories. The 
realizations that sputtering occurs primarily in the first 
three layers of the target and that focusons do not, toa 
great extent, affect the sputtering yield were made in this 
manner. This thesis will use the MI program QDYN86 to 
‘simulate the problem outlined below. 

While a computer simulation can provide insight into the 


sputtering process, it is gwst that; a  smmularuoene an 
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tieterotandine OF the phvysreal sttuation must be Attained and 
balanced by physics and mathematics before a simulation can 
be successfully achieved. Some fundamentals, such as the 
interatomic potential Bnet 1OnS wa cabeomm stall snot willy 
understood. In this respect, computer simulations cannot and 
will not supercede experimentation. They can however, 
provide support and explanation of observed phenomenan = and, 
mere importantly, ae used as a cutting edge of science in 


the probing of new ideas and concepts. 


D. APPLICATIONS 

Sputtering simulations are not limited to monocrystalline 
targets. Indeed, much of the effort currently being expended 
lies in the arena of alloys and chemically reacted surface 
bombardment. Chemical reaction of the surfaces involves the 
introduction ot an additional, different atom onto the 


surface of a crystalline target. This additional atom is 


called an adatom. The locations of these adatoms, privsr 75 
bombardment, is usually unknown, however, tor various 
crystalline structures, anticipated locations may be 


described. Both tungsten and molybdenum are BCC lattices 
Sreuehures and the four-fold position (Pig 1) is a possible 
candidate, depending on the amount of coverage desired. In 
that location, the adatom would sit in a position between 
the four face atoms and a center atom such that either the 


forces acting on it, such as lattice and surface binding 
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energies, equalize, or it could be in a position such that 
the distance to it’s nearest neighbors would be equal (Fig 
1). Whether the adatom is more likly to be located above or 
below the surface plane and which of the above situations is 
applicable is an unanswered question for most systems. 

The realization that important real world applications 
can be addressed in Simulations has spurred their 
application to surface science investigations. Corrosion, 
the formation of oxides on the surface of a substance, can 
be studied. Possible solutions to its destructive ettects 
such as adatom placement to either prohibit or slow down the 
introduction of the oxygen atom into the lattice can be 


Simulated. Conversly, positioning may allow the adatom to 


act as a catalyst, enhancing a desired reaction. 
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Il. OBJECTIVES 


Winters (Ret 20) has conducted experiments investigating 
sputtering cross sections of both tungsten and molybdenum 
targets which have been reacted with nitrogen adatoms, when 
bombarded by various ions. These experiments showed that the 
nitrogen sputtering yield tended to increase as the atomic 
weight of the target atoms increased. Thus, assuming a (001) 
face on the target, the sputtering yield of the nitrogen 
adatoms was reasoned to be dominated by the mass ratio of 
the substrate to the adatoms. Meyerhoff (Ref 21), using a 
Simulation model, determined that the distance between the 
absorbate and the substrate atoms was more significant than 
the mass ratio. 

Both Winters and Meyerhoff assumed that the location of 
the nitrogen adatoms was in the four-fold position (Fig 1). 
Meyerhoff conducted his simulation assuming a nitrogen 
location slightly above the suface plane, resulting in a 
position otf equal distance between the adatom and all five 
nearest-neighbors. Earlier work by Calvenna and Schmidt (Ref 
22) and Adams and Germer (Ref 23) also assumed the four-fold 
position however, there was little regard given to target 
atom displacement due to the introduction of the nitrogen 
adatoms. Griffiths et. al. (Ref 24) incorporated this idea 
in a series of experiments using low energy electron 


diffraction (LEED). Using the observed patterns for various 


1S) 


coverages and primary beam energies, they were able toa show 
that, at a fractional coverage ( defined as the ratio, for a 
specific face of the target, of the number ot adatoms to the 
number otf substrate atoms.) of 0.5 and a surface temperature 
of 1000°K, the nitrogen adatoms were unstable and were 
rapidly absorbed into the bulk. If the coverage was reduced 
to approximatly 0.4, they found that temperatures a te 
ITOC0SK di denen Beenie in any further absorption. A coverage 
of 0.4 however, required the formation of "islands" of 
substrate with nitrogen adatoms on the target (Fig <2). The 
target used in this thesis can be interpreted as one oft 
these "islands". To investigate this type of structure, 
various laser diffraction gratings were constructed and 
tested by Griffiths et al. A grating with a random number of 
nitrogen adatoms, but with an average island size of 4X4, 
produced a pattern similar to that observed when LEED beams 
were used on an actual target although the intensities of 
various parts of the pattern were not similar. Harrison (Ref 
25) suggested that the nitrogen adatoms could be located 
below the surface close to an interstitial site. As stated 
earlier, Meyerhoff conducted his simulation assuming the 
equilibrium position of the nitrogen adatom above the 
surface plane of the substrate. This thesis will study 
equilibrium points both above and below the substrate 
surface plane. 


The primary objective of this thesis is to study 
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sputtering from BCC crystalline targets of molybdenum = and 
tungsten which have had nitrogen adatoms placed in different 
equilibrium positions, above and below the substrate surface 
plane for both molybdenum and tungsten targets. The 


incident ions will be argon of various energies. 


ae 


III. Simulation Model Development 


A. QDYN8&5 

The simulation program used for this thesis research 
was QDYN85. This program is an MI Simulation using a time 
step approach. Initial inputs of crystalline structure, 
atomic masses, potential function parameters, adatom 
locations, ion eee eee and impact point are required. Using 
Hamilton’s equations of motion, the program initiates and 
develops the subsequent cascade through the crystal until 
the energy of the atoms is lowered to the point where no 
futher ejections can occur. 

After the initial ion impact, POS1 ti one velocity 
directions and energies of affected target atoms and adatoms 
are computed every timestep. To compute these values ftor 
every atom, every timestep, the actual program running time 
(cost) would be excessive. Therefore, forces are not 
computed on a particular target atom until it is actually 
hit by a moving atom. Additionally, the program maintains a 
listing which is updated periodically, of each atom’s 
nearest neighbors. 

For this particular simulation, the timestep is variable 
and is determined by a specified distance Gdiviced bya 
highest atomic velocity in the target. The distance choosen 
for this calculation is 0.1 lattice units (LU). A lattice 


unit is defined as one half the lattice parameter of any 
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GUbMEMeCL Sealine Structure mine lattice unit was Created to 
facilitate calculations concerning crystal structure bv 
eliminating the need to input a large number of new 
variables when there is a change in the elements being 
studied. Lo determine the appropriate timestep more 
accurately, the rough “distance over velocity” calculation 
1s modified by limiting Pacuere which take into account the 
previous velocity pane ensuring a smoother transition in 
time. 

Harrison and Jakas (Ref 26) present a more detailed 
description and analysis of the MI code of QDYN85, as well 
aS a comparison of its results with those obtained ‘by 


TReiSer, a binary collision simulation. 


B. POTENTIAL FUNCTIONS 

Of critical importance to the sputtering process is the 
modeling of how the atoms of the crystal interact with each 
other. To date, no single potential function has been 
developed that fully and accuratly reflects every sputtering 
experiment. While newer potential functions may be developed 
or the parameters changed to fit the observed results, the 
basic mechanisms remain stable. Four basic potential 
functions are used here. 

1.Born-Mayer 

This is a strictly repulsive function which is 


considered adequate for intermediate separations. It has the 


io 


form: 
VirJ= a exp 40- or e 


2.Moliere 


This is also a repulsive function but of the form: 


Vir)=(Z1 42e4/r)(0.35exp(-0_ 3277 a) -Umose se be 


+0 .lexp(-6.0r/a)], 
where "a" is the Firsov length and is calculated as tollows 
a=0:8853anr GA 2472 ee le oe 


where ab is the Bohr radius. Because the Moliere function 
was originally intended as an approximation Ge) the 
Thomas-Fermi function,it is accepted that the value obtained 
for "a" may be modified. This initiated a new terminology, 
the "modified Moliere" potential function. 
3.Morse 

The Morse is a potential function which is both 

attractive and repulSive depending on the separation 


distance ‘r". It has the following form 


V(r)=Deexpl=2&(r-re )})} = 2beexpl 2 seam 
where De is the well depth, re is the equilibrium separation 
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of an atomic pair and alpha is a parameter controlling the 
Shape of the well 
4.Morse-Moliere 

Tacs function iS a composite potential function 
consisting of sections of both the Morse and the Moliere 
functions. The two are joined together at a specific section 
by a cubic spline. The intersection is determined by varying 
the alpha and "a" parameters of the Morse and Moliere 
functions respectivly to obtain an intersection where the 
Slopes are as close together as possible. The result is a 
smooth joining with the repulsive wall governing the 
collision dynamics and an attractive well which governs the 
sputtering. 

Obviously, values obtained by these functions are 
dependent on the elements used for the target, adatom and 
the incident ion. Table 1 gives the pertinent data for the 
elements used in this simulation. The value “e is the 
lattice parameter, B is the binding energy of nitrogen to 


the substrate and C is the coheSive energy. 


Table 1: Natural Data for Elements Used. 


At.Mass Structure ao (A) B(eV) CleV) 


Mo gO 4 BCU eee | N/A 6.820 
W SSeS D BCG Oreo 5 N/A 8.900 
Ar 6 ..0 N/A N/A N/A N/A 
N 14.0 N/A N/A One N/A 


5. Substrate-Substrate Function 


The potential function used for the crystalline 


at 


structure is a composite . Morse-Moliscre Pune Lola The 
development of the potential parameters for this interaction 
will be essentially the same for the substrate-adatom, 
substrate-ion and adatom-ion interactions described below. 
As stated earlier, the alpha and the ‘a" values were 
choosen to match the respective potential graph slopes. The 
value of De was choosen so that the energy from the atoms of 
the lattice, oneneoG at the center of the lattice, equaled 
the cohesive energy listed in table 1. The value of Re is 
the nearest neighbor separation distance. The values of Ra 
and Rv,the spline boundaries, were determined by the final 
Sraphs of the potentials which, in turn, were dictated by 


the alpha and "a" values. There is also a parameter Re which 


is used in the program. It’s function i185 #topelamt tie 
distance, in lattice units, over which, the force and 
potential calculations are made. This value is 1.71 #42xLU 


except for the substrate-substrate interaction in which case 
Re is 2.40 LU. 
6. Adatom-Substrate Function 

The parameters for this function are determined in 
the same manner as the substrate-substrate function with the 
exception of the Re and De values. In this instance, the 
value Re is the measured separation distance vertically 
between the cubic center and the adatom and De is the 
binding energy of the nitrogen atoms to the surface, 


6.5eV/atom. 


Ze, 


t. lTon-Substrate/Adatom Functions 


In addition to these adatom-substrate potential 
functions, others must be considered. Since nitrogen has 
been introduced into the target, the possibility exists 
that, at some time, two nitrogen atoms may become neighbors 


or even collide. The potential function used to describe 
this evolution is a pure Morse function with the following 


experimentally determined values: 


Sas 2. 700 Ge-# , 
Re= 1.098 8, 
“DOES TI eV), 


During the simulation of a cascade, a single incident ion 
is injected onto the target. Since it is unlikly that the 
argon ion will react with the nitrogen, molybdenum or the 
tungsten, these potential functions are asSigned as 
“modified” Moliere’s. 

The overall initial governing factors in the 
determination of the above values in this’ and the two 
previous sections were the desire to obtain a smooth force 
curve and for the minimum point of the well depth to occur 
when the adatom was below the surface plane. Tables 2 and 3 
show the comparison between the values assuming an adatom 
location of 0.05 LU below the crystalline plane, and 


Meyerhoff’s. Table 4 shows the values assuming various above 


aa 


plane positions but incorperating the interatomic forces. 


Table 2: POTENTIAL FUNCTION PARAMETERS. 
(PREVIOUS) 
De (eV) Re (A) &(8-1) MOD a Ra (LU) Rb (LU) 
N-Mo 1, G2) 149 50Nee aca 910 970 
N-W 1.300 2.360" 1 e5@ mane 1.1407 dan 
Mo-Mo .997 2.800 1.500 0.0 790 830 
W-W 1.335 2.849 1.200 0.0 1.830 1.900 
Table 3: POTENTIAL FUNCTION PARAMETERS (NEW) 
(BELOW PLANE) 
De (eV) Re (A) cy(A-1) MOD a Ra(LU) Rb(LU) 
N-Mo 2.970 1.524 23600 - 80g Se 560 
N-W 3.292 1.524 2.850 -.800 °#®.650 670 
Mo-Mo 9o7—demeCOnna 5 om One 790 830 
W-W 1.335 2.894 1.200 0.0 1.120. 3. 
TABLE 4: CRYSTAL POTENTIAL FUNCTION PARAMETERS 
ABOVE PLANE LOCATION 
LOCATION De Re o< MOD a Re Rp 
(eV) (A) (&-!) (nu)  “Caup 
Mo-N .380 1.600 1.950 2.150 0.0 720 750 
Mo-N .245 1.997 1.738 2.143 -.8 500 520 
W-N .487 1.712 1.960 2.320 g 830 850 
Figures 3-10 show the potential and force graphs 
for the new values of the crystal potential function as 


compared to those used previously by Meyerhoff. A review of 


the force graphs using Meyerhoff’s data, shows a possible 


discontinuity point in the area where the splining between 


the 


Morse and Moliere potentials occurs. It should be noted 


tuhat, at the time of the research, the effects the 


interatomic forces generated would have on the final results 


was not recognized. The differences are not significant. 


Figures 11 through 14 #=show the plots for the old and new 


Mo-Mo and W-W forces. Figure 15 shows the graph of the 
nitrogen-nitrogen potential function and figure 15 shows the 
Eraphs of the ion-adatom and the ion-sSubstrate potential 


funetions. 


8.Vacuum Phase Potential Parameters 

The probability exists that interactions above the 
surface of the eee will occur. The potentials used for 
these interactions are all pure Morse. There is a paucity of 
data for the Mo-Mo, W-W, Mo-N and W-N systems and, as_ such, 
approximations had to be made. The values required are De, 
Re , “omega” and “omega-x" where the later two are functions 
of the vibrational and rotational energy of the system. The 
value @ is also required and 1s calculated, after 


determining De, as follows: 


K 


(1.36116x10-3 )(omega)(m/De )1/2 


m = reduced mass of the system 


Mo data could be found for the Mo-Mo, W-W, Mo-N or the W-N 
systems. Therefore, aS a compromise, using reference 27, 
systems were found that were as close to the investigated 
system as possible in reduced mass and in location on the 
periodic table. For Mo-Mo, Ag-Br was used and ~VO was 
substituted for Mo-N. Similar substitutions were used for 


the tungsten system. The N-N system was available. Table 5 


Zo 


summarizes the values used. 


TABLE 5: VACUUM PHASE POTENTIAL PARAMETERS. 


De Re omega omega-x 
Mo-Mo 1.700 1.730) 2.393024 ce .679 
N-N 5.760 3.359 1.03908 2350557/0 ‘143225 
Mo-N 6.100 1,965. 1.52457 0Z0meo0 4.700 
W-W 3.100 174650" °2- 470 ~3037502 . 960 
W-N 4.900 2.159 1/38 2967008 4.850 


C. TARGET and IMPACT AREAS 

The size of the target is determined by the energy of the 
incident ion. If the selected target size is insufficient, 
it is possible that, during the course of a cascade, atoms 
which have escaped from the sides or the bottom of the 
target would have caused additional sputtering to occur if 
the target had been larger. Such an event is calle a 
“failure of containment”. Although a large target is 
desirable to maximize the sputtering yield, there exists a 
practcal upper bound. For a given incident energy, 
increasing the size of the target beyond some point would 
waste computer resources. Theretore, a balance must be 
struck between the two requirements. Using a copper target 
bombarded with argon, it was found that a target size of 
23X8X23 for an incident energy of 2 keV would provide almost 
complete containment. 

For this size target, the introduction of 60 nitrogen 


adatoms, each in the four-fold position, would result ina 
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Peace b2ona weoverarecmonumee!)/; so Pieureei, shows the (001) face 
for such a target. Although nitrogen is the adatom used in 
this thesis, other elements may be used. This could result 
in the adatoms occupying a different location on the target. 
These locations may be the two-fold (equidistant between two 
nearest-neighbor substrate atoms) or the A-top (directly 
above a substrate atom) positions (figure 18). 

A total of 600 impact points on the surface were examined 
for each incident ion energy value. These impact points, as 
measured on the x and z axis, were varied about ae central, 
pre-selected Dolne on the target. There were two 
mre-selected points, (11.000, 11.000) ana@ (13.000, Loo O } 
(Fig 17). The first location does not have an occupying 
adatom while the second does. This method results ina large 
number of independent runs for a specific ion energy and 


therefore, more accuracy in the results. 
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IV. RESULTS AND DISCUSSIONS 


A. ANALYZING AND PLOTTING PROGRAMS 

ANMOL is the analyzing program currently being used fosr 
sputtering research at the Naval Postgraduate School. The 
results obtained from QDYN&5 are used as input with the 
ANMOL program formatting and tabulating the results. The 
output can be varied aS to the depth and detail desired fer 
& given input. 

The inputs used for the ANMOL oroeban are also those used 
tor the program ANPLOT which will provide graphical vice 
tabular results. In addition to ANPLOT, the graphics program 


EASYPLOT was used to generate all force and potential plots. 


Bo Hill /NO alt SCHINAR TO 

. As discussed previously, each Simulation run was 
conducted with a total of 300 impact points referenced to a 
pre-selected point on the target surtace. For the 
Simulation, only half of the available positions contain a 
nitrogen adatom. Therefore, there is an equal probability of 
the incident ion hitting or missing an adatom. To 
incorperate this difterence, runs, at the same ion energy, 
were conducted for both instances where the ion hit and did 
not hit a nitrogen adatom with the results being averaged 
together. This results in essentially a 600 point run for 


each energy simulated. The primary difference of the two 
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runs is that there are approximatly two to four times the 


number ot sputtered nitrogen atoms for the hit run as there 


are tor the no-hit run, depending on the ‘incident. ion 
energy. 
All table listings ot yield -and cross sections 


incorperate this hit/ no-hit scenario. 


C. NITROGEN RESULTS» 
1.Simulated cross sections and yields 

The value of the sputtering yield of nitrogen can be 
expressed as the product of the concentration of the 
nitrogen adatoms in atoms per square centimeter (@n ) and 
wre Sputtering Cross section in sdWare centimeters (On ). 
The results are in units of nitrogen atoms per ion. Figure 
Moeeeshows the experimental results of Winters for the 
sputtering cross section of nitrogen versus the ion energy 
for a variety of substrates. 

The simulated concentration of the nitrogen, assuming a 
coverage of 0.4, was calculated to be 5.0EF14 atoms/cm2. 
Table 6 shows the results of the nitrogen cross sections for 
both molybdenum and tungsten substrates with the adatoms 
located 0.05 LU beneath the surface plane. 

Since the adatom/substrate potential functions change 
for different adatom positions, different parameters must be 
calculated for each assumed position of the adatom to ensure 


equilibrium. New potential function parameters were computed 


for various adatom positions above the plane for both the 
molybdenum and tungsten systems. Table 7 shows the results 
£ or the selected positions (in parentheses below the 
substrate type) above the plane for each system. The unit 


for this distance is the LU. 


TABLE 6: SIMULATION NITROGEN CROSS SECTIONS 
CROSS ,SECTION X10-1!15(cm2) (BELOW PLANE) 


ION ENERGY YIELD CROSS SECTION 
(KEV ) Mo(001) W(001) Mo(001) W(001) 
1 .33 N/A . 66 N/A 
2 mae 29 (oe eof 
3 er N/A 46 N/A 


TABLE 7: SIMULATION NITROGEN CROSS SECTIONS 
CROSS SECTION X10-15 (cm2) (ABOVE PLANE) 


ION ENERGY YIELD CROSS SECTION 
(KEV ) Mo(001) W(O01) Mo(001) W(O01) 
(Pago... Jae) (.487) (.245 .380) {346m 

i N/A oa -2Ou N/A 1.14 lode 

Zs 30M oe . 450 .6 OS 2 I 

J N/A oe N/A N/A 1.10 N/A 


2. Nitrogen results comparison to experimental results 
Figure 20 shows a comparison for the molybdenum 
substrate with adatom positions of 0.05 LU below the surface 
for the low point and 0.38 LU above the surface for the high 
point. 
When the cross section and yield simulation results for 
the molybdenum target with the below plane adatom location 


are compared to the experimental values obtained by Winters, 


a0 


a strong negatively sloped curve is evident rather than a 
armen poasttlvs elope Also, tie  Similation “cross section 
values are below the experimentally observed values. This 
indicates that the adatom probably is not located beneath 
the surtace of the target. The observed cross section for 
the below plane nitrogen from tungsten bombarded at 2 keV 
was obviously low, suggesting that the adatom was not 
located below the aetnee Tor that system either. This 
observation negated the need for further runs on tungsten at 
this adatom moet tiOn. 

Figure 20 also compares the computed molybdenum results 
Ho x the two previously mentioned above plane adatom 
locations to those obtained by Winters. With respect to the 
0.38 LU position, the curves are much closer in shape, and 
the new points tor molybdenum are now slightly above the 
experimental values. This seems to indicate that the assumed 
position of 0.38 LU above the plane was too high a location 
for the adatoms. 

The equidistant positioning of the adatom, which was a 
change of 0.135 down LU from the previously examined point, 
resulted in a decrease in the cross section of 40%. This 
large change in yield for a small change in distance 
supports the distance dependency postulated by Meyerhoff. 
These results also indicate that the position of the 
nitrogen adatom, for a molybdenum target, is slightly above 


the eqidistant position. Figure 21 and shows a comparison of 


on 


Uae experimental values obtained for tungsten substrate and 
those abtained from the simulation for the "above plane’ 
adatom location of 0.487 LU. These appear to be in fairly 
close agreement. However, little can be definitivly stated 
Since only two points were evaluated for the assumed 
POs fikon. 
3. Simulated nitrogen energy distribution 

Figures 22 PO oueh 20 show the energy distribution of the 
sputtered nitrogen atoms for an incident energy of 2 keV. In 
aaa case, the spectrum was fairly broad. In both extreme 
cases, where the adatom was located at it’s closest and 
furthest position from the surface, for both tungsten and 
molybdenum, a spectrum with a peak of approximaly 5 eV is 
evident; therefore, the choice of the substrate, when there 
is a large mass difference between the substrate and the 
adatom, appears to be of little consequence to the energy 
distribution of the ejected nitrogen. Additionally, the 
location of the adatom has little effect on it’s sputtered 
energy. These results conflict with those obtained by 
Meyerhoff. 
4.Simulated nitrogen ejection angles 

Figure 26 is a spot pattern presentation of the polar and 
azmuthal angle dependence of the sputtered nitrogen adatoms. 
The azmuthal angle of the spot is read normally while the 
polar angle of a spot is plotted as the numerical value of 


the tangent of the polar angle. Therefore, spots further out 


a2 


Ke DAK 7 K 
; v A J es 


angles of ejection as 


with larger 


adatoms 


correspond to 


any , 


at 


Little, 


measured from the surface normal. 


pattern 


in 


seen from the graph. Table 8 lists, 


be 


can 


development 


typical polar and azmuthal angular results 


numerical format, 


for sputtered nitrogen atoms of various energies. 


Nitrogen atom souttering angles 


Table 8: 
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D. MOLYBDENUM AND TUNGSTEN RESULTS 
When the incident ion impacts a target with a chemisorbed 
gas, target atoms as well as the adatoms are sputtered. 
Clean surface targets of both molybdenum and tungsten 
were bombarded in the simulation with 2 keV argon ions. The 
number of sputtered target atoms were then compared to the 
number of sputtered target atoms when the target had been 


reacted with nitrogen. Table 9 lists these values. 


TABLE 9: SUBSTRATE ATOM EJECTIONS 


CLEAN POALOMeEOS I? TON 
HIGH LOW 

Mo 907 677 686 ~ 
W ol 70 So) 


A comparison of the number of sputtered target atoms for 
the low and high nitrogen positions of both substrates was 
done. The comparison shows that although the nitrogen acts 
aS an energy absorption mechanism, it’s location is of 
little consequence to the ejection mechanism of the 
substrate material. 

1. Simulation target yield 

Figures 27 and Z8 show the target atom yield per impact 
point for a 2 keV incident ion. The bottom left corner of 
the diagram is the selected. reference point for the ion 
impact. Throughout the graph, the darker the point, the 
higher the yield for that particular point. The lighter 
areas of the graphs are indicative of the channelling effect 


for the BCC structure. 
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a. Simulation Target Yield Comparison to Experiment 

In his experiments, Winters also calculates the 
ratios of the cross sections of the chemisorbed gas to the 
substrate material. Table 10 lists the results of the 
Simulation as compared to the experiment. Winters used 
polycrystalline tungsten and molybdenum targets in his 
experiments while the calculations of this thesis used clean 
target (no adatoms) cross sections. The simulation tungsten 
and molybdenum concentrations were calculated to be 1.19E15 


and 1.2E15 respectivly. 


TABLE 10: CROSS SECTION RATIO’S WITH NITROGEN IN EQUIDISTAN 


POSITION 
EXPERIMENT SIMULATION 
JON ENERGY CN(W) ON(MO) ~N(W) o~N(MO) 
O-W( POLY) C”MO( POLY) O W(CLEAN) © MO(CLEAN) 
c . wae .68 .41 24 


Since the simulation clean targets were bombarded with zd keV 
ions only, comparisons must be restricted to that energy. 
The simulation cross section value for nitrogen from the 
molybdenum target was calculated assuming the 0.245 position 
for the adatom. Table 10 shows a smaller value than the 
espeimentar Sara although the nitrogen cross section ratio 
for experiment and simulation are both 1.33. These results, 
along with the previously mentioned Significant 
yield/distance relationship indicate that the assumed 


Positions are close to the actual adatom locations. 


“Iona .On smoaanoadbilty of ejection 

Figures 29 and 30 are representitive of the probability 
Simo JCCETCa ee tor Gach of the surface atoms oan this case on 
the (001) face. These figures indicate how well the cascade 
sequence was contained by the target. The graphs show that 
the size of the target was generally acceptable for 
molybdenum but was’ slightly small for tungsten. This does 
not afmeet the epee results for either the substrate 
or the adatom significantly. 
3. Simulation substrate ejection spot patterns 

Figures 31 and 32 #£show typical sputtered target atom 
spot patterns. As with the nitrogen pote, these sive 
graphical representations of the angles of ejection of the 
target atoms. The crystalline nature of the targets is 
readily apparent aia these fraphs | because of the 
concentration of sputtered atoms along the 45°9,135°,225°0, 
and 315° arms. For a BCC structure, these directions 
correspond to the nearest neighbor directions 
4.Simulation ejected substrate energy distribution 

Figures 33 through 38 show typical energy distribution 
histograms of the sputtered target atoms. There were no 
apparent differences between the specific substrate results. 
The molybdenum results show a maximum at approximatly 5eV 
while the maximum for tungsten is slightly higher. These 
values conflict with the theoretical values which should be 


one halt the cohesive energy; specifically, 3.41 eV for 


36 


molybdenum and 4745 eV tor tungsten. When plotted 
logarithimicly, both elements approximatly show the expected 
E-2 dependency at higher energies (Fig.39 through 43). Table 
11 shows the slope relationship numerically with the values 


having been obtained from the figures mentioned. 


TABLE 11: SUBSTRATE EJECTED ATOM ENERGY DISTRIBUTION 


SLOPES 
ION ENERGY SUPSTRATE ENERGY 
(KEV) DISTRIBUTION SLOPE 
3.0 MO lee 
1.0 W -1.93 
2.0 W -1.94 
270 W (clean) -1.99 
Ze MO (clean) = ei 
TS 


E. DIMBR RESUL 
| When a compound structure is bombarded, there may be an 
expectation that groups of two or more atoms would be 
Sputtered off together. The results of this simulation do 
not show this to be the case. In many instances however, 
after ejection, two atoms, either substrate - substrate or 
Substrate - adatom combine over the target. Individual 
ejection times of combining atoms showed differences of as 
much as 200 fsec, indicating that combination occured at a 
relativly large distance from the target. However, there 
were some instances when the times of ejection were 
separated by as little as 10 fsec for two atoms whose 
pre-ejection locations were within 1 LU thus indicating 


dimer formation almost immediatly after ejection. 


ot 


VY. CONCLUSIONS AND RECOMMENDATIONS 


The results of the simulations, when compared Lo 
experimental results, indicate + hist , how nitrogen 
chemisorbded on the (001) face of molybdenum and tungsten 
targets, the location of the adatom is above the surface 
plane for both targets. 

Nitrogen yield Sanna ons between the simulation and 
Winters experimental results indicate that the adatom 
position for the molybdenum target is above the equdistant 
position of 0.245 LU and below the 0.380 LU position. 
Furthermore, the simulation yield comparisons between the 
two positions show that there is a very strong dependency of 
the yield on the adatom distance from the target surtace 
plane. 

The comparison of the simulation nitrogen yield to the 
experimental values for the tungsten target indicate that 
Mieeeadabtom poSition iS near the equidistant position of 
0.487 LU. However, since only two ion energy values were 
analyzed for this target, further investigation of the 
tungsten target should be conducted at different ion 

energies to complete the argon simulation results. 
| hMaaachonavivy, s  turther Simulation investigation into both 
molybdenum and tungsten should be conducted using different 
bombarding ions in the same energy range as those used in 


this thesis. These additional simulations should be analyzed 


38 


in conjunction with the results obtained in this thesis to 
isolate the exact adatom locations for the targets in 


question. 
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omall circle represents nitrogen atoms. 


Fig. 1. Location of nitrogen atoms on target substrates 
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The proposed contracted-domain structure 

for the 0.4 monolayer on a Mo(001) or a 

W(001) surface: Large hatched circles, top- 
layer substrate atoms, small filled circles 

N atoms. 

a) Plan view, illustrating domain and boundary 
structure. 

b) Cross section through the domain along AA: 


Fig. 2. Island formation on a target surtace 
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Tungsten-Nitrogen force function used previously 
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for tungsten target 
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Fig. 23. Ejected Nitrogen energy distribution from a 
molybdenum target assuming argon energy 2 keV 
and adatom location -.05 LU 
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Fig. 24. Ejected nitrogen energy distribution from 
molybdenum target assuming argon energy of 2 keV 
and adatom position .387 LU 
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Fig. 25. Ejected nitrogen atom energy distribution from a 
simulated tungsten target assuming argon energy 
of 2 keV and adatom position .487 LU 
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Fig. 29. Simulation atom ejection probability from a 
PCC (001) molybdenum target assuming a 2keV 
argon ion 


2.000 W(001) +60N/A<001> (23X8X23) TUNGSTEN (NH) 
ATOM EJECTION PROBABILITY 














BCC (001) 
0.096O.0066O0 
€-0-0:0-0:0-6-6-0-6"° 
OHH CRC On OER 
tet tet Bod Boe 
Ct SR © 
eee Se 
7 geese. 
oe: Goon 
COG Pog 
x0 ( “) ©) 

OH) ox see (COM (CO) 

OxCrEUG 

oe greet xc 

Soe 
0:6: 0-6:0-6:0-650-6 

Ge AEE EE CHCHER 

SOOOUSOOOOOOUL 


Fig. 30. Simulation atom ejection probability from a 
BCC (001) tungsten target assuming an argon 
ion energy of 2 keV 
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Fig. 32. Simulation tungsten ejection spot pattern assuming 
2 keV argon ion 
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Fig. 35. Simulation ejected atom energy distribution from 
tungsten target with nitrogen adatoms assuming a 
1 keV argon ion 
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Fig. 37. Simulated ejected atom energy distribution from 


clean molybdenum target assuming a 2 keV argon 
ion 
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Fig. 38. Simulated ejected atom energy distribution from 
clean tungsten target assuming a 2 keV argon ion 
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Fig. 39. Log-Log ejected atom energy distribution from a 
molybdenum target with nitrogen adatoms assuming 
3 keV argon ion 
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Fig. 40. Log-Log ejected atom energy distribution from 
tungsten target with nitrogen adatoms assuming 
1 keV argon ion 
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Fig. 41. Log-Log ejected atom energy distribution from 


tungsten target with nitrogen adatoms assuming 
2 keV argon ion 
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Fig. 42. Log-Log ejected atom energy distribution from 
clean tungsten assuming 2keV argon ion 
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Fig. 43. Log-Log ejected atom energy distribution from 
clean molybdenum assuming 2 keV argon ion 
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